Composition-induced structural transitions in mixed rare-gas clusters 



F. Calvo 

Laboratoire de Physique Quantique, IRSAMC, Universite Paul Sabatier, 
118 Route de Narbonne, F31062 Toulouse Cedex, France 

E. Yurtsever 

Kog University, Rumelifenenyolu, Sariyer, Istanbul 34450, Turkey 

The low-energy structures of mixed Ar-Xe and Kr-Xe Lennard- Jones clusters are investigated 
using a newly developed parallel Monte Carlo minimization algorithm with specific exchange moves 
between particles or trajectories. Tests on the 13- and 19- atom clusters show a significant im- 
provement over the conventional basin-hopping method, the average search length being reduced by 
more than one order of magnitude. The method is applied to the more difficult case of the 38-atom 
cluster, for which the homogeneous clusters have a truncated octahedral shape. It is found that 
alloys of dissimilar elements (Ar-Xe) favor polytetrahedral geometries over octahedra due to the 
reduced strain penalty. Conversely, octahedra are even more stable in Kr-Xe alloys than in Kr3g 
or Xe38, and they show a core-surface phase separation behavior. These trends are indeed also 
observed and further analysed on the 55-atom cluster. Finally, we correlate the relative stability of 
cubic structures in these clusters to the glassforming character of the bulk mixtures. 



I. INTRODUCTION 

Clusters of heterogeneous materials show a much richer 
behavior than their homogeneous counterparts. In many 
bulk compounds, doping can significantly affect some 
global property, and alloying is a common way to tai- 
lor a completely new kind of material. At the mesoscale 
level, size is another complicating factor, giving rise to 
further changes with respect to the macroscopic object. 
To a large extent, most expectations of nanotcchnology 
have been put into the electronic and catalytic properties 
of small atomic clusters. Therefore, it should not be sur- 
prising that numerous theoretical studies of mixed clus- 
ters were devoted to bimetallic clusters. In particular, 
there has been a significant amount of work at the level 
of sophisticated electronic structure calculationSfii2i£ but 
these were often limited to small sizes due to the nu- 
merical effort involved. On a different scale of chemical 
complexity, many studies have been carried out using 
explicit, empirical force fj e lds^*Si2i2i2iiSiiiiiS in order to 
investigate the segregation properties of these clusters. 

There are several driving forces toward mixing or seg- 
regation in binary systems: 

(i) the difference in atomic sizes; 

(ii) the difference in surface energies; 

(iii) minimization of the overall strain; 

(iv) the number of interactions between unlike atoms. 

These factors can often compete with each other. For in- 
stance, minimizing surface energies does usually not in- 
crease with the number of interactions between different 
atoms. Also, even though this is not our prime interest 
here, it should be noted that kinetic factors can be crucial 
in this problem^ 

In particular, Vach and coworkers have found from ex- 
periments and simulations of mixed rare-gas clusters that 



some anomalous enrichment effects could be observed 
due to the growth by pick-up of these systems^ Very 
recently, radial segregation and layering have been ob- 
served in large Ar/Xe clusters formed in an adiabatic 
expansion by Tchaplyguine et al. 15 using photoelectron 
spectroscopy measurements. These data have also been 
theoretically interpreted by Amar and Smaby^ 

Fortunately, mixed rare-gas systems can be quite safely 
described using simple pairwise potentials such as the 
Lennard-Jones (LJ) potential. More accurate potentials 
are of course also available, even though we will have 
no need for them in the present, mostly methodological 
work. Hence they are much more convenient to study in 
a broad size range, not only for their structure but also 
their dynamics or thermodynamics. It is known from 
previous studies that the topography of the potential en- 
ergy surfaces of homogeneous LJ clusters can be very 
peculiar, as for the sizes 38 or 75 The multiple-funnel 
structure of these energy landscapes makes it especially 
hard to locate the most stable structures (global min- 
ima) or to simulate the finite-temperature behavior of 
these clusters in an ergodic way. The effects of mixing 
different rare-gas atoms on cluster structure and ther- 
modynamics have been studied for the specific size 13 
by Frantz on the examples of Ar-Kr mixtures^ as well 
as Ne-Ar mixtures^ Fanourgakis et al. have also in- 
vestigated these latter compounds^ A systematic work 
of Ar-Xe mixed clusters of 13 and 19 atoms has been 
carried out by Munro and coworkers, 21 including some 
global optimization and Monte Carlo simulations. Mixed 
clusters involving lighter species such as H 2 and D 2 have 
been investigated using path-integral Monte Carlo simu- 
lations (PIMC) by Chakravarty^ More recently, Sabo, 
Doll and Freeman reported a rather complete study of 
the energy landscape a 23 ! 24 a nd melting phase change 2 ^ in 
mixed Ar-Ne clusters. In this work quantum derealiza- 
tion and the effects of impurities on cluster properties 
were also accounted for using PIMC techniques. 
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The main conclusion of these studies is that atomic 
heterogeneity can be responsible for a drastic increase in 
complexity of the energy landscapes of rare-gas clusters. 
This complexity is manifested by numerous new low- lying 
minima in competitive funnels, characterized by the same 
overall geometrical arrangement but different permuta- 
tions of unlike atoms. Following Jellinek and Krissinel, 4 
we will refer to such isomers as "homotops." The pres- 
ence of several homotops on a given energy landscape of- 
ten induces solid-solid transitions, which can be detected 
by some feature in the heat capacityji2*i2i22i2iiSi even 
though they can be washed out by quantum effects^ 
As shown by Munro et a/.f2i the various funnels corre- 
sponding to different homotops of a same geometry are 
separated by significant energy barriers. This explains 
the difficulty or even failure of simulation methods to 
achieve ergodic sampling of these systems, albeit small 
A similar situation is found in Lennard- Jones polymers, 27 
where a large number of isomers are based on the same 
geometrical arrangement, differing only in the path link- 
ing the monomers. 

Beyond the actual rare-gases, binary Lennard-Jones 
compounds have been investigated in both the cluster 
and bulk regimes. Clarke and coworkers looked at phase 
separation of small particles with equal compositions. 28 
Based on Monte Carlo simulations, they sketched a 
phase diagram in the general structure of liquid clus- 
ters. Bulk binary Lennard-Jones systems have been seen 
to provide relatively simple numerical models for glass 
formationi 29 i 30 i 31 i 32 i 33 i 34 Most often, the LJ interactions 
in such studies have been tuned in a non-additive way in 
order to hinder crystallisation. In another related work, 
Lee and coworkers 3 ^ have investigated the role of atomic 
size ratio in binary and ternary metallic alloys. 

Interestingly, severaly links between the physics and 
chemistry of clusters and those of supercooled liquids 
and glasses have been established since the pioneering 
work by Franks The initial suggestion that the local or- 
der in simple liquids is not crystalline but icosahedral 36 
(more generally polytetrahedral) has since been veri- 
fied experimentally 37 and theoretically. 29138 From the 
clusters viewpoint, the favored finite-size structures of 
good model glassformers have been shown by Doye and 
coworkers to be polytetrahedral^ 9 . 

The 38-atom homogeneous Lennard-Jones cluster is 
known to show some glassy properties, especially slow 
relaxation to the ground stated due to the competi- 
tion between two stable funnels on the energy landscape, 
corresponding to truncated octahedral and icosahedral 
shapes, respectively. Due to entropic effectsj^Sdi a solid- 
solid transition occurs between the two funnels, at tem- 
peratures lower than the melting point. The crystal- like 
configuration of this cluster makes it a good candidate to 
further investigate the relationship between cluster struc- 
ture and criteria for glassification. 

Because homogeneous L J38 constitutes a relatively dif- 
ficult task for global optimization algorithms, binary 
clusters of the same size can be expected to be much 



worse. In this paper, we propose a simple but efficient 
way to deal with the multiple new minima introduced by 
unlike atoms within a general Monte Carlo global min- 
imization scheme. This algorithm will then be applied 
to the 38- and 55-atom cases, in mixtures of Xe with ei- 
ther Ar or Kr atoms. In the next section, we present the 
method and test it on the simple cases of the 13- and 
19-atom clusters. In Sec. lIIII we give our results obtained 
at sizes 38 and 55 and we correlate them to the different 
glassforming abilities of the bulk mixtures. We finally 
conclude in Sec. IIVI 



II. METHODS 

Global optimization of cluster structured is cur- 
rently best achieved using either genetic algorithms^ 3 - or 
the Monte Carlo+minimization method^ also known 
as basin-hopping (BH).~ The case of homogeneous 
Lennard-Jones clusters is among the most documented 
of cluster physics, and an up-to-date table of putative 
global minima can be found in Ref. l46t Even though 
it can never be guaranteed that global minimization has 
been successful, it is likely that all important structural 
forms of LJ clusters have been found up to more than 100 
atoms. These include icosahedral, truncated octahedral, 
decahedral as well as tetrahedral arrangements. 

Compared to homogeneous clusters, the available data 
on heterogeneous systems is rather scarce. Besides the 
specific works by Frantz on the 13-atom Ne-Ar and Kr- 
Ar clustersjiSii 9 . Munro et al. used a parallel version of 
the BH scheme, similar to the replica-exchange Monte 
Carlo method^ where several trajectories are run simul- 
taneously at various temperatures. 21 Although these au- 
thors looked at moderately large clusters, they reported 
significant difficulties to locate global minima at specific 
compositions, as in XeioAr3 or XeiaArg, for instanced 



A. Optimization algorithm 

A natural problem occuring using the basin-hopping 
method is that many of the low-lying minima arc ex- 
pected to be related to each other via particle exchange. 
Such a process only occurs via large deformations of 
the remaining cluster, hence it is quite unprobable. As 
in condensed matter physics ^ILMLSiL allowing exchange 
moves between particles as a possible Monte Carlo step 
may result in notably faster convergence, provided that 
the interactions are not too dissimilar. Actually, op- 
timization of mixed clusters on the lattice formed by 
the homogeneous system has already been studied by 
Robertson and coworkers*^ Here we do not wish to re- 
strict to such situations. 

In the framework of global optimization methods, the 
local minimization stage removes the possible energetic 
penalty associated to replacing a small atom by a big- 
ger one. We can thus expect some increased efficiency 
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of the algorithm in case of multiple homotops. Now we 
convert the extra numerical cost of running parallel tra- 
jectories at various temperatures into running them at 
various compositions, at the same fixed temperature T 
for all compositions. For a X p Y„_ p compound, each of 
the n trajectories is then labelled with the number p of 
X atoms, running from to n. Exchange moves between 
adjacent trajectories (from p to p + 1) thus need to in- 
corporate the transmutation of two atoms (one for each 
configuration) into the other atom type to preserve com- 
position. As in most Monte Carlo processes, the proba- 
bility of attempting such moves must be set in advance 
as a parameter. 

The global optimization algorithm can thus be summa- 
rized into its main steps. Keeping the above notations 
for atom types, and denoting Tl^ the configuration at 
step i of trajectory p, we start the optimization process 
using fully random configurations, but locally optimized. 

1. With probability P cx , it is decided whether an ex- 
change between adjacent trajectories will be at- 
tempted or not. If so, then the two trajectories 
involved in the exchange are determined randomly. 

2. For each composition p not concerned by any ex- 
change, a new configuration is generated from 

(p) 

R^ using either several particle exchanges or large 
atomic moves. The probability to select particle ex- 
changes is denoted P S wap , and the number of simul- 
taneous exchanges is allowed to fluctuate randomly 
between 1 and iV^J. If atomic moves are selected, 
then each atom is displaced randomly around its 
previous location in the three directions by a ran- 
dom amount of maximum magnitude h^ p \ In both 
cases, Rj^x is obtained after local minimization. 

3. In case of an exchange between adjacent trajecto- 
ries, the two configurations and R^ p+1 ^ cor- 
responding to these trajectories are then swapped, 
one X atom of R^ being transmuted into Y, and 
one Y atom of R| p+1 ' ) being transmuted into X. 

Again, the configurations R.[?_i and R^ 1 "* are ob- 
tained after local minimization. 

4. Each new configuration is accepted with the usual 
Metropolis acceptance probability at temperature 
T. 

The algorithm has two main parameters, namely P cx 
and P SW ap- The maximum number of particle exchange 
moves, -/V s ™|p, was set to 4 in this study. We expect 
that better results could be obtained by adjusting this 
parameter appropriately, probably taking higher values 
for larger clusters or for compostitions close to 50%. The 
amplitude of atomic displacement, h^ p \ is set to half the 
equilibrium distance in the X2 dimer for p = 0, half the 
equilibrium distance in the Y2 dimer for p = n, and is 
interpolated linearly between these two values for < 



p < n. In the present work, the exchange probabilities 
were taken as P cx = 0.5 and -P S wa P = 0.9, hence allowing 
a rather large probability of sampling among homotops 
of a same structure. 



B. Benchmark calculations 

Low-energy structures for mixtures of xenon with ei- 
ther argon or krypton atoms have been first investigated 
for the sizes n = 13 and n — 19, as there are quantita- 
tive global optimization data available for Ar-Xe clusters 
from the Jordan group. 2 - We have adjusted the LJ values 
used by Leitner et alJH to reproduce the clusters ener- 
gies found by Munro and coworkersiSi. With respect to 
argon, the present data for a and e are thus <7KrKr = 
1.12403, axexe = 1.206, a Kr xo - 1.16397, <T ArX c = 1-074, 
e&Kr = 1.373534, exoXo = 1.852, e KrX e = 1.59914, and 
£ArXc = 1-48. Global optimization of Ar-Xe and Kr-Xe 
clusters was performed using the parallel algorithm pre- 
viously described, simultaneously for all compositions, 
for a maximum number of 10000 minimization steps per 
trajectory, and at T = 0. Ten independent runs were 
carried out to estimate an average search length for each 
composition. All global minima reported by Munro et 
al. were always found within the number of MC steps 
allowed. 

The results for Ar„Xei3_„ an d Kr„Xei3_„ clusters are 
given in Tabled The average search length is generally 
higher for compositions close to 50%, for which the num- 
ber of homotops is maximum for a given isomer, regard- 
less of symmetry. The statistics presently obtained for 
Ar-Xe clusters show that the average search is between 
10 and 1000 times faster than using conventional paral- 
lel basin-hoppingiSi Kr-Xe clusters roughly exhibit the 
same level of difficulty, but we do not see any strong ev- 
idence for particularly severe cases: AraXeio even seems 
to be one of the easiest. 

Similarly, the results obtained for Ar^Xeig-n clusters 
show a significant improvement over fixed-composition 
basin-hopping.— They are given in Table [D] along with 
the corresponding data for Kr„Xeig_„ clusters. This 
time, the algorithm is about 1-100 times faster depending 
on n, the average search length being still longer for equal 
compositions. For both the 13- and 19-atom clusters, all 
global minima are homotops of either the single or dou- 
ble icosahedron. This situation is particularly suited for 
our algorithm, especially the exchange moves. 

Initially, the configurations at all compositions are ran- 
dom. The chances to locate the proper structure (with- 
out any consideration of the homotops) increase linearly 
with the number of trajectories. As soon as the right 
structure is found, the algorithm naturally optimizes 
atom types to find the most stable homotop, hence the 
global minimum. But it can also communicate the struc- 
ture to the adjacent trajectories, until all compositions 
only need to sample among the permutational homotops. 

When the interactions are not too dissimilar (as in Kr- 
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TABLE I: Global optimization results for Ar„Xei3_ n and Kr n Xei3_„ clusters. The search length is the average over 10 
independent runs of the number of Monte Carlo steps needed to find the global minimum. Energies are given in LJ units for 
argon. 



Ar n Xei3_ro 


Global minimum 


Average 


Kr n Xei3_ n 


Global minimum 


Average 


cluster 


energy 


search length 


cluster 


energy 


search length 


Xei3 


-82.093 


3.2 


Xei3 


-82.093 


3.0 


ArXei2 


-78.698 


7.9 


KrXe i2 


-81.014 


9.6 


Ar 2 Xen 


-76.274 


9.6 


Kr2Xen 


-79.263 


4.3 


Ar-^Xei n 


-74.015 


5.8 


Kr^Xei n 


-77.550 


5.7 


Ar4Xeg 


-71.597 


8.6 


Kr4Xeg 


-75.869 


26.1 


ArsXeg 


-69.017 


14.0 


KrsXes 


-74.186 


25.8 


Ar6Xe7 


-66.584 


37.2 


Kr 6 Xe 7 


-72.498 


26.4 


Ar 7 Xe 6 


-63.791 


19.4 


Kr 7 Xe 6 


-70.844 


45.2 


Ar 8 Xe 5 


-60.733 


13.9 


Kr 8 Xe 5 


-69.141 


18.3 


ArgXe4 


-57.851 


22.7 


KrgXe4 


-67.473 


4.7 


Ari Xe 3 


-54.594 


12.0 


Kri Xe 3 


-65.802 


11.5 


ArnXe2 


-51.122 


7.5 


KrnXea 


-64.128 


4.1 


AriaXe 


-47.698 


4.1 


KriaXe 


-62.490 


2.4 


An 3 


-44.327 


2.7 


Kri 3 


-60.884 


2.3 



TABLE II: Global optimization results for Ar n Xei9_ n and Kr n Xei9_ n clusters. The search length is the average over 10 
independent runs of the number of Monte Carlo steps needed to find the global minimum. Energies are given in LJ units for 
argon. 



Ar„Xei9_„ 


Global minimum 


Average 


Kr„Xeig_„ 


Global minimum 


Average 


cluster 


energy 


search length 


cluster 


energy search length 


Xei 9 


-134.566 


72.4 


Xeig 


-134.566 


70.7 


ArXeig 


-131.819 


64.3 


KrXeig 


-133.651 


94.0 


Ar2Xei7 


-129.116 


80.3 


Kr2Xei 7 


-132.701 


109.8 


Ar 3 Xei6 


-126.547 


85.2 


Kr 3 Xei 6 


-130.088 


84.3 


Ar 4 Xei5 


-123.764 


238.2 


Kr 4 Xei 5 


-129.067 


167.2 


Ar 5 Xei4 


-120.786 


196.6 


KrsXei4 


-127.284 


175.4 


Ar 6 Xei 3 


-118.284 


221.2 


Kr 6 Xei3 


-125.498 


265.9 


Ar 7 Xei2 


-115.681 


391.8 


Kr 7 Xei 2 


-123.709 


334.6 


ArgXen 


-113.075 


387.9 


KrgXen 


-121.951 


319.5 


Ar 9 Xeio 


-110.242 


264.2 


KrgXeio 


-120.115 


287.1 


AnoXeg 


-107.531 


295.8 


KnoXeg 


-118.304 


243.6 


AmXeg 


-104.576 


193.8 


KmXeg 


-116.521 


187.4 


Ari2Xe 7 


-101.811 


235.5 


Kri 2 Xe 7 


-114.736 


214.3 


An 3 Xe 6 


-98.110 


158.5 


Kn 3 Xe 6 


-112.947 


201.3 


Ari 4 Xe 5 


-94.396 


247.3 


Kn 4 Xe5 


-111.189 


188.8 


Ari 5 Xe 4 


-90.438 


121.3 


Kri 5 Xe 4 


-108.863 


176.5 


Ari 6 Xe 3 


-86.328 


131.2 


Kri 6 Xe 3 


-106.609 


115.1 


Ari 7 Xe2 


-81.907 


97.3 


Kri 7 Xe 2 


-104.332 


10.2 


ArigXe 


-77.298 


86.8 


KngXe 


-102.036 


98.1 


Ar 19 


-72.660 


62.2 


Kng 


-99.801 


69.8 



Xe clusters), it is likely that the mixed clusters share the 
same isomer as the global minimum of the homogeneous 
cluster, which justifies the lattice approach of Robertson 
et idmiL The problem is then reduced to locating the most 
stable homotops. By setting P S wap to one and starting 



all trajectories from this minimum, the algorithm can 
be even more successful, and we estimated the average 
search length to be further reduced by a factor about 3 
with respect to the values given in Table [HI However, 
when the interactions differ significantly among atoms 
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types, or when the energy landscape of the homogeneous 
cluster does not display a single steep funnel, it becomes 
much harder to make a guess about structure in these 
binary clusters. 

4000 




10 

Xe atoms 



FIG. 1: Mean first passage time of the parallel optimiza- 
tion algorithm to locate the global minimum structure of 
Ari9_ n Xe n clusters versus n. The average is performed over 
10 independent runs. The results are shown at zero temper- 
ature, with or without exchange moves (E) between trajecto- 
ries, with or without swap moves between atom types (S). The 
results at T = 1 with both kinds of moves are also displayed. 

Fig. ^ shows the mean first passage time needed to 
locate the global minima of Arig_ n Xe„ using the algo- 
rithm under different conditions. Disabling swap moves 
between atom types or exchange moves between adja- 
cent trajectories usually attenuates the efficiency. Em- 
ploying a rather high temperature is even worse, because 
the cluster may easily leave its optimal lattice. This con- 
trasts with optimizing homogeneous clusters, where using 
a nonzero temperature helps the system to escape from a 
funnel. '' However, if the energy gap between homotops 
of the same lattice increases and gets close to the gap be- 
tween different lattices, we expect the zero temperature 
method not to be the best. But in such cases, even the 
notion of a lattice should be questionned. 



III. STRUCTURAL TRANSITIONS 

In this section we focus on two larger sizes, for which 
no global optimization result is available. The LJ38 
cluster is characterized by its archetypal two-funnel en- 
ergy landscaped The high free-energy barrier separating 
these two funnels and the higher entropy of the less sta- 
ble minima of the icosahedral funnel make it particularly 
hard to locate the truncated octahedral lowest-energy 
minimum using unbiased global optimization algorithms. 



Hence it is not surprising that this peculiar structure was 
first found by construction. 56-57 



Composition-induced transitions in the 38-atom 
clusters 



We have attempted to locate global minima for binary 
Ar-Xe and Kr-Xe clusters of size 38, using the parallel 
basin-hopping algorithm previously described. Because 
of the huge number of homotops at this size, and most 
importantly because of the structural competition be- 
tween icosahedra and truncated octahedra, we cannot be 
fully confident that the global optimization was success- 
ful. Therefore, the energies reported in Table ITTT1 for Ar- 
Xe clusters should be taken with caution, as they could 
probably be bettered. The same data for 38-atom Kr-Xe 
clusters is also reported in Table Ivl 

Specifically to this cluster size, all minima found dur- 
ing the optimization process were categorized as either 
icosahedral or cubic-like, depending on the energy of the 
corresponding homogeneous isomer found by quenching. 
In cases where the cubic isomer was not found among the 
isomers, we performed additional optimizations starting 
from this structure, setting P S wap to one. This mainly 
occured for Ar-Xe clusters. Eventually, two series of 
minima were obtained for each of the icosahedral and 
octahedral funnels. We did not find any decahedral iso- 
mer that could compete with these structural types, even 
though some evidence for stabilizing decahedra by doping 
was reported in Ref. |S3 



LU 

< 




20 

Xe atoms 



FIG. 2: Energy difference AE = Ef cc — E[ co between the most 
stable cubic and icosahedral isomers of Kr38_ n Xe n (empty 
circles) and Ar38-nXe n (full squares). AE is given in reduced 
Lennard- Jones unit of Argon (approximately 120 K). 

We have represented in Fig. [3 the relative energy dif- 
ferences AE = Ef cc — Eico between the most stable cubic 
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TABLE III: Global optimization results for Ar38-nXe n . The energies are given in LJ units for argon, and the symmetry and 
mixing ratios defined by Eq. Q are reported. 
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isomers and the most stable icosahedral isomers, as they 
were obtained from our optimization scheme, for both the 
Ar-Xe and Kr-Xe mixtures. Besides some strong varia- 
tions sometimes seen from one composition to the next, 
and which can be attributed to usual finite-size effects, 
general trends can be clearly observed. 

First, Kr38_„Xe„ clusters are always most stable in the 
cubic shape. Actually, changing the composition most 
often further stabilizes truncated octahedra, and only 
rarely enhances the stability of icosahedra, which occurs 
for ?i > 29 and n = 21. Conversely Ar 38 _ n Xe n clus- 
ters are preferentially found icosahedral, exceptions being 
n > 35 and n = 0. This is an example of a composition- 
induced structural transition between the two funnels of 
the energy landscape. 

From a computational point of view, it should be noted 
that the optimization algorithm was able to locate the 
truncated octahedral minima for Kr-Xe clusters by it- 
self, starting from disordered minima, and that the extra 
runs starting from this structure only produced slightly 
more stable homotops. This is another illustration of the 
efficiency of the present parallel optimization method. 

The general degree of disorder is higher in icosahedral 
structures than in the cubic-like isomer. Hence it is more 
difficult to put up the latter geometry with very unlike 
interactions, as in Ar-Xe clusters. Cubic homotops of 
argon with xenon are rather distorted, but the strain is 
much lower with krypton instead of argon. Examples of 
global minima obtained at compositions n = 9, 19, and 
29 are represented in Fig. [3] In Kr-Xe compounds, a 



progressive core-surface phase separation is seen with Kr 
atoms outside, in agreement with energetic arguments: 
atoms with the larger e prefer to occupy interior sites. 
In icosahedral clusters, the strain increases at such sites, 
especially in polytetrahedral systems. Icosahedral Kr-Xe 
clusters also prefer to have Xe atoms at the center, but 
the increased strain is too high a penalty, which explains 
that cubic structures are favored over icosahedra. 

In general, no complete phase separation is found in 
Ar-Xe clusters, even though Ar atoms seem to fit best at 
the centre of the cluster. In both cases, surface energies 
thus play an important role. Mixing in these clusters can 
be estimated using radial distribution functions^* Here 
we use the same index as Jellinek and Krissinel^ namely 
the overall mixing ratio 7 defined asm 

j(X p Y n -p) = 
Ex p y„- p - E Xp (X p X n - p ) ~ E Yn _ p (YpY n -p) 

where Ex p Y n _ p is the binding energy of cluster X p Y n - Pl 
Ex (XpX n ^p) the binding energy of subcluster X p in the 
homogeneous cluster X p X n - p at the same atomic config- 
uration as XpY n - p , and a similar definition for the last 
term of Eq. As seen from Table ITTT1 the mixing ratio 
increases notably in Ar-Xe clusters, up to more than 60% 
for some compositions. Kr-Xe clusters, despite exhibit- 
ing some core-surface segregation, show similar variations 
of the mixing ratio with composition, with only slightly 
smaller values of 7. Therefore the mixing ratio, as de- 
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Ar g Xe 2g (Ci) Ar 19 Xe 19 (C x ) Ar 29 Xe 9 (C 2 „) 




Kr 9 Xe 29 (C„) Kr 19 Xe 19 (C,) Kr 29 Xe 9 (C 4 „) 



FIG. 3: Putative global minima found for several Ar38- n Xe„ 
and Kr38- n Xe n clusters. Argon, krypton, and xenon atoms 
are represented as black, gray, and white balls, respectively. 
The point groups are indicated. 

fined in Eq. QJ, is a rather ambiguous parameter for 
quantifying the extent of mixing in this small cluster. 

The optimal structure of an homogeneous cluster de- 
scribed with a pairwise potential results from a compe- 
tition between maximizing the number of nearest neigh- 
bors and minimizing the strain energy, or penalty in- 
duced by distorting these bonds, 5 ' 7 Binary Lennard-Jones 
systems exhibit several extra complications due to the 
various ways of rearranging atom types in a given struc- 
ture. In these systems, the strain varies notably among 
the homotops, especially in clusters made of very unlike 
atoms. However, our results indicate that the same gen- 
eral rules hold for homogeneous and heterogeneous sys- 
tems. In ArigXei9, the strain is rather high, but the num- 
ber of contacts is also high. In KrigXei9, both the strain 
and the number of nearest neighbors are much smaller. 

To investigate the role of heterogeneity on the strain, 
we have computed the various contributions to the re- 
duced strain energies in Ar38_„Xe n clusters. The strain 
energies are defined for Ar-Ar, Xe-Xe, and Ar-Xe inter- 
actions as follows^ 7 - 

^strain T/LJ , AT nn 

^Ar-Ar — ^Ar- Ar + ly Ar- Ar £ Ar-Ar , 

restrain T/LJ , AT nn 

-^Ar-Xc — ^Ar-Xc + JV Ar-Xc £ ArXe, 

Tristram taLJ i AT^n- ^~ 

-^Xe-Xo - l-Xc-Xc + iV Xe-Xe £ Xo-Xc- 

In these equations, V^_ Y i s the (negative) total binding 
energy between atoms X and Y, is the number 

of X-Y nearest neighbors, and £x-y is the Lennard- 
Jones well depth corresponding to the interaction be- 
tween X and Y atoms. Reduced strain energies are then 
defined as e strain = E 3tTain /N nn e, in order to account for 
the different magnitudes of the interactions among atom 



types. According to these definitions, the strain ener- 
gies are always positive quantities. The strain energies 
in Ar38_„Xe„ clusters are represented versus composition 
in Fig. They give us some insight about the possible 
ways of reducing strain. 



0.15 




Xe atoms 

FIG. 4: Reduced strain energies for alike and unlike interac- 
tions in Ar38-nXe„ clusters, versus composition n. 



The pattern exhibited by the reduced strain versus 
composition shows different behaviors for clusters having 
mostly argon or xenon atoms. For n < 20, most strain 
is carried by interactions between alike atoms. This 
case is illustrated by ArgXe29 m Fig- 01 where a kind 
of core/surface phase separation occurs. Here surface 
energies are also important, but the situation is rather 
different from mixed cubic Kr-Xe clusters. Because hav- 
ing the xenon atoms at the inner sites of the icosahedral 
structure would maximise the strain of these atoms, it 
is much more favorable to have the smaller atoms inside 
and the xenon atoms outside. The cubic Kr-Xe struc- 
tures, on the other hand, are not especially strained, and 
having the smaller atoms inside would lead to an ener- 
getic penalty. 

When the number of Ar atoms increases above about 
19 in the 38-atom cluster, interactions between unlike 
atoms are significantly more strained. The case of 
Ar2gXeg depicted in Figgis perticularly informative: Xe 
atoms are located scarcely among the icosahedral cluster, 
and relieve the structure from too much strain at the ex- 
pense of only few Xe-Xe interactions. In this case, cluster 
structure is driven by the number of unlike interactions. 

It is also worth noting that a few compositions are 
especially weakly strained; this occurs when the global 
minimum is octahedral, but also in the range 19 < n < 
24. For these latter clusters, the core/surface segregation 
and the number of unlike interactions are both optimal. 
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Polytetrahedral transitions in the 55-atom 
Ar— Xe clusters 



The cubic to icosahedral transition seen above actu- 
ally favors polyicosahedral (or anti-Mackay) structures. 
The strain reduction produced by size disparity in 38- 
atom Ar-Xe clusters helps in stabilizing these kinds of 
structures, which are otherwise replaced by multilayer (or 
Mackay) geometries in the homogeneous clusters. Since 
most L J clusters under the size of 38 atoms are most sta- 
ble as polytetrahedral we do not expect that changing 
composition will affect them to a large extent. As a no- 
table exception, the 6-atom homogeneous LJ cluster is 
more stable in its octahedral isomer. The lowest energy 
geometries of mixed Ar-Xe clusters containing 6 atoms, 
represented in Fig. [5J show polytetrahedral transitions 
for two compositions, namely Ar4Xe2 and Ar3Xe3. This 
behavior mimics somewhat what was observed for the 
larger 38-atom cluster, only at a smaller scale. In par- 
ticular, and as in Fig. |3 polytetrahedral arrangements 
are seen to be more convenient for Xenon compositions 
under 50%. 

Possible polytetrahedral structures of mixed Ar-Xe 
clusters have been investigated for the size 55, whose 
most stable isomer is well known as a perfect double 
layer (Mackay) icosahedron for the homogeneous system. 
The global optimization results are summarized in Ta- 
ble IIVI for all compositions. For this we also conducted 
complementary calculations on the 2-layer lattice. Each 
putative global minimum was labelled either as Mackay 
icosahedron or, when the lattice structure does not ex- 
actly match the multilayer icosahedron, as polytetrahe- 
dral. Most compositions become increasingly polytetra- 
hedral as the ratio of Xenon atoms increases, even though 
the polytetrahedral character may often be only local. 

Two examples of lowest energy structures are repre- 
sented in Fig. H3 corresponding to n = 15 and n — 40. 
The obvious deviations of the geometry of the former 
from the Mackay icosahedron and the various occupa- 
tion sites of the heavy atoms for both structures illus- 
trate again that there is no simple rule that determine 
the most stable minima when the atomic sizes do signif- 
icantly differ from each other. 



C. Temperature-induced transitions 

We now go back to the 38-atom clusters of Kr and Xe 
atoms, for which the global minimum was always found to 
be a truncated octahedron. The extremely large number 
of isomers (including homotops) in the energy landscape 
of binary Lennard- Jones clusters, added to the expected 
presence of significant energy barriers between icosahe- 
dral and cubic isomers^ prevent finite-temperature sim- 
ulations from being conducted in a reliably ergodic way 
with the presently available tools. For example, the par- 
ticle exchange moves used to accelerate convergence of 
sampling among homotops will likely have very low ac- 



ceptance probabilities in MC simulations at low tempera- 
tures, especially for Ar-Xe clusters. Therefore, even with 
powerful methods such as parallel tempering or multi- 
canonical Monte Carlo, reaching convergence in 38-atom 
LJ clusters does not seem currently feasible to us. 

As an alternative, we have chosen to investigate 
solid-solid transitions by means of the superposition 
approaches! For a given cluster, databases of minima in 
each of the icosahedral (ICO) and truncated octahedral 
(FCC) funnels were constructed using the optimization 
algorithm. For each composition and each of the two fun- 
nels, no more than 2000 distinct minima were considered. 
The classical partition function of the Y38_ p Xe p cluster 
(Y=Ar or Kr) restricted to funnel A=FCC or ICO is ap- 
proximated by a harmonic superposition over all minima 
of the databases, which belong to this funnel: 60 



Qa(P) 



(2) 



where f3 — l/kgT is the inverse temperature, i>i the ge- 
ometric mean vibrational frequency, = 2p\(n — 
with hi the order of the point group of minimum i and 
n = 38. We do not consider quantum effects here, al- 
though they may be important at low temperatures jSi 
since derealization or zero-point effects are not expected 
to be significant for rare gases as heavy as krypton or 
xenon. 

Within the harmonic superposition approximation, a 
solid-solid transition occurs when Qfcc = QiCO~ This 
latter equation is solved numerically in (3 or T, its so- 
lution is denoted T ss . In cases where icosahedra are 
energetically more stable than octahedra, a solid-solid 
transition can occur if some cubic structures are entrop- 
ically favored, which requires lower vibrational frequen- 
cies and/or lower symmetries. We did not find such sit- 
uations in our samples of Ar-Xe clusters, therefore we 
restrict to Kr-Xe clusters in the following. 

Similar to transitions between funnels, transitions be- 
tween homotops will happen if their partition functions 
are equal. The huge number of homotops gives rise to as 
many values for the corresponding temperatures, and we 
define the homotop transition temperature Th such that 



T h = mm{T^\T^ } > 0}, 



(3) 



where TP' is the transition temperature between the 
global minimum (homotop 0) and its homotop j. 

Equating the harmonic partition functions for these 
two isomers leads to the expression of T^'tr^ 



Ej — Eq 



(3n — 6) In v§ /i/j + In rij / no 



(4) 



Since all homotops are characterized by different vibra- 
tional and symmetry properties, the transition temper- 
atures are not ordered exactly as the energy differ- 
ences Ej—Eo. This reflects that solid-solid transitions in- 
volve crossover in free energy rather than binding energy. 
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FIG. 5: Lowest energy stru 




Ari 5 Xe 40 Ar 40 Xei 5 



FIG. 6: Lowest energy structures of the ArisXe4o and 
Ar4oXei5 clusters. Both structures have Civ symmetry. 



The above equation also shows that can take nega- 
tive values if homotop j has a higher symmetry and/or 
a higher vibrational frequency than the ground state. In 
this case the global minimum is always the free energy 
minimum, and no solid-solid transition occurs, hence the 
form of Eq. © . 

Finally, a third temperature has a strong consequence 
on cluster structure, namely the melting temperature. 
Its estimation from either simulations or superpositions 
approximations is already quite difficult for the homo- 
geneous LJ 38 cluster ; 40 ! 41 and we did not attempt to 
compute it for binary clusters. However, the previous 
study by FrantaiS has shown that the melting point 
in mixed, 13-atom Ar-Kr clusters varies quite regularly 
(approximately quadratically) with composition. As a 
simple rule, we will assume that the melting point of 
K r 38-nXe„, T me it(n), lies inside some range between the 
approximate melting points of KT38 and Xe38, respec- 
tively. From the results obtained by Doye and Wales 4 ^ 
and the Monte Carlo data of Ref. |4l] for the LJ38 clus- 
ter, we get T mc it(0) ~ 0.234 and T melt (38) ~ 0.315 in 
reduced LJ units of argon. This provides rough limits to 
the actual melting points of Kr-Xe clusters, for the price 
of neglecting finite-size effects. 
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FIG. 7: Solid-solid transition temperature (full circles) be- 
tween the octahedral and icosahedral funnels, estimated from 
a harmonic superposition approximation, versus composition 
in Kr38_„Xe n clusters. Also shown is the lowest transition 
temperature for a permutation between octahedral homotops. 
The horizontal dotted lines mark upper and lower limits for 
the estimated melting points. The temperatures are given in 
LJ units for argon (120 K). 

The transition temperatures are represented in Fig. [7| 
for all compositions in Kr38_„Xe„ clusters. We no- 
tice first that the structural transition temperature T ss 
varies quite regularly with composition in both the ranges 
n < 19 and n > 21, and that it shows strong size effects 
between these limits. Several situations are predicted to 
occur depending on the relative values of T ss , Th and 

Tmclt ■ 

In most cases, T me i t < T S s- That melting takes place at 
temperatures lower than the cubic/icosahedral transition 
simply nullifies the transition between structural types. 
However, this extra stability of the octahedral funnel may 
have a consequence on the melting point itself, which is 
likely to increase. Still, this situation implies that simu- 
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TABLE IV: Global optimization results for Ar55_„Xe n - The energies are given in LJ units for argon, and the symmetry and 
mixing ratios defined by Eq. Q are reported. The structural types (Mackay Icosahedron MI or polytetrahedral PT) are also 
given. 



n 


Mixing 


ratio Energy 


Point group Type 


n 


Mixing 


ratio Energy 


Point gr 


oup Type 








-279.248 


1 n 


MI 


28 


0.53 


-417.785 


Cx 


PT 


1 


0.05 


-284.276 




MI 


29 


0.53 


-421.225 


C 3 


PT 


2 


0.10 


-289.313 


c 

w s 


MI 


30 


0.51 


-424.566 


Cx 


PT 


3 


0.15 


-294.360 


c 


MI 


31 


0.49 


-427.888 


Cx 


PT 


4 


0.28 


-302.344 




PT 


32 


0.47 


-431.565 


Cx 


PT 


5 


0.34 


-310.780 


c 


PT 


33 


0.46 


-435.680 


Cx 


PT 


6 


0.38 


-316.331 


c 


PT 


34 


0.45 


-439.565 


Cx 


PT 


7 


0.41 


-321.770 




PT 


35 


0.42 


-443.335 


Cx 


PT 


8 


0.43 


-327.241 


c 


PT 


36 


0.41 


-446.892 


Cx 


PT 


9 


0.44 


-332.356 


d 


PT 


37 


0.41 


-449.965 


Cx 


PT 


10 


0.47 


-337.497 




PT 


38 


0.36 


-454.356 


C 3 


MI 


11 


0.50 


-342.655 


c 


PT 


39 


0.35 


-458.203 


C s 


MI 


12 


0.55 


-347.608 


d 

^- 1 


PT 


40 


0.34 


-462.564 




MI 


13 


0.57 


-352.520 


d 


PT 


41 


0.30 


-466.709 


Ch 


MI 


14 


0.56 


-357.586 


c 


PT 


42 


0.28 


-472.191 


h 


MI 


15 


0.58 


-362.483 




PT 


43 


0.26 


-475.967 




MI 


16 


0.58 


-367.231 


d 

^ 1 


PT 


44 


0.24 


-479.739 


D 5d 


MI 


17 


0.59 


-372.016 


C- 


PT 


45 


0.22 


-483.495 




MI 


IS 


0.59 


-376.740 


d 


PT 


46 


0.20 


-487.219 


C'2v 


MI 


19 


0.59 


-381.501 




PT 


47 


0.18 


-490.910 


C 3 


MI 


20 


0.60 


-386.007 


Cs 


PT 


48 


0.16 


-494.585 


Ci 


MI 


21 


0.60 


-390.135 


Ci 


PT 


49 


0.14 


-498.231 


Cs 


MI 


22 


0.59 


-395.065 


C x 


PT 


50 


0.12 


-501.860 


C2v 


MI 


23 


0.59 


-399.213 


C x 


PT 


51 


0.10 


-505.467 


Civ 


MI 


24 


0.60 


-403.205 


Cx 


PT 


52 


0.08 


-509.052 


D 5d 


MI 


25 


0.58 


-406.983 


Cx 


PT 


53 


0.06 


-512.616 




MI 


26 


0.56 


-410.686 


Cx 


PT 


54 


0.04 


-516.170 


h 


MI 


27 


0.56 


-414.427 


C a 


PT 


55 





-517.168 


h 


MI 



lations will more easily reach convergence. 

However, there are notable exceptions for this behav- 
ior, at n < 4, n = 21 and n > 34. In these clusters, 
heterogeneity is not sufficient for the thermodynamical 
behavior of the cluster to deviate too much from those of 
the homogeneous system. 

The transition between homotops usually occurs prior 
to melting. Thermal equilibrium within the cubic funnels 
thus involves several homotops (and "restricted" solid- 
solid transitions), and the corresponding thermodynam- 
ical state could be probably simulated using specifically 
designed exchange moves between outer particles within 
a Monte Carlo scheme. 

A few clusters melt before exhibiting any transition 
between homotops. This occurs for instance at n = 11, 13 
or 15. For these sizes the structural transition also occurs 
at temperatures higher than the estimated melting point. 
These cases should pose less problems to conventional 
simulations than the homogeneous cluster. 



D. Glassy behavior 

The previous results have shown that finite-size Ar- 
Xe compounds show a preferential polytetrahedral order, 
even for very low doping rates, over octahedral order. On 
the other hand, Kr-Xe clusters at the same sizes further 
favor cubic order. Since polytetrahedral order is known 
to be present in liquids and, more generally, in disor- 
dered structural glasses, it seems natural to correlate the 
behavior observed in these clusters to the dynamics of 
the corresponding bulk materials^ 

We have simulated the cooling of 108-atom binary 
rare-gas liquids, using a simple Metropolis Monte Carlo 
scheme under constant volume and temperature. Ini- 
tially the atoms are placed randomly into a cubic box of 
side L, and periodic boundaries are treated in the min- 
imum image convention. The LJ interactions were not 
truncated, and the simulations consisted of 100 stages of 
10 5 MC cycles each, linearly spaced in temperature. 

Three compositions have been selected, following our 
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knowledge of the cluster structure. For each composition, 
different length sizes L and different temperature ranges 
Proin i Tmax] were chosen in order to cover both sides of 
the melting point. In the first mixture, 24 xenon atoms 
and 84 argon atoms are simulated with L = 4.8815 LJ 
units of argon, with 0.1 < T < 1. In the second mix- 
ture, 24 xenon atoms are added to 84 krypton atoms at 
L = 5.487 and 0.15 < T < 1.5. The third mixture con- 
sists of 9 argon atoms and 99 xenon atoms at L = 5.887 
and 0.2 < T < 2. Even though we did not attempt to 
locate the most stable crystalline forms for these mix- 
tures, our searches close to the face-centred cubic mor- 
phology showed that the most stable configurations for 
these mixtures always had some cubic order. It is likely 
that the actual ground states for such systems are indeed 
crystalline^ 

The average root mean square fluctuation of the bond 
distances, also known as the Lindemann index S, uni- 
versally characterizes the thermodynamical state of the 
condensed system as either solid or liquid, depending 
on its value being lower or higher than about 0.15. To 
quantify the extent of crystalline order, we have used the 
bond order parameter Q4 introduced by Steinhardt and 
coworkers^ The two parameters S and Q4 allow us to 
follow in Monte Carlo time the cooling processes for all 
materials in a simultaneous way, independently of ther- 
modynamical characteristics such as the melting temper- 
ature. 

0.20 1 : 1 . , . , . 1 
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FIG. 8: Correlation between the Lindemann parameter S and 
the order parameter Q4 along cooling simulations of Ar-Xe 
and Kr-Xe bulk mixtures. 

The correlation between S and Q4 for ten cooling sim- 
ulations of each of the three bulk binary compounds is 
represented in Fig. [8] In all cases, the Lindemann param- 
eter covers the whole range 0.01 < S < 0.18, indicating 
that the melting point was indeed crossed. However, the 
three compounds display very contrasted cooling behav- 
iors. 



In the (Arg4Xe24) system, 5 regularly decreases but 
Q4 always remain below 0.05. Therefore crystallization 
never takes place, and the final state obtained by quench- 
ing is significantly higher in energy than some crystalline 
forms; this is typical of glass formation. 

In (Ki'84Xe24), all simulations show some rather sharp 
transition from a (high 5, low Q4) state to a (low 5, 
high Q4) state as 5 crosses about 0.1. The tempera- 
tures where crystallisation occurs may vary somewhat 
among the cooling runs, in the same way as they are ex- 
pected to depend on the cooling rate. Lastly, the case of 
(ArgXegg) is intermediate: while most simulations end 
up in a nearly fully crystalline phase (Q4 ~ 0.15), a few 
of them show a limited degree of cubic ordering in the 
solid phase, Q4 having values close to 0.07. 

These results very closely reflect our previous data on 
binary, 38-atom clusters of the same materials. In terms 
of composition, the first mixture corresponds to Ar3oXes, 
which clearly favors icosahedral shapes over truncated 
octahedra. The second mixture reminds of Kr3oXe§, for 
which the cubic structure is even more stable than in 
the homogeneous cluster. The third mixture should be 
compared to Ar3Xe35, which favors icosahedra only mod- 
erately. 

This correlation found here between cluster structure 
and the glassforming ability of the bulk material con- 
firms previous analyses on the icosahedral local order in 
liquids and glasses , 29 i 38 as well as the recent conclusions 
obtained by Doye et alS^ that clusters of good glassform- 
ers indeed show a polytetrahedral order. 



IV. CONCLUSION 

As far as structural and dynamical properties are con- 
cerned, binary compounds show a significantly richer 
complexity with respect to homogeneous clusters. The 
work reported in the present paper was intended to 
achieve several goals. First, a parallel global optimiza- 
tion algorithm was designed to locate the most stable 
structures of mixed rare-gas clusters, beyond the lattice 
approximations of Robertson and co-workers^i Based 
on the basin-hopping or Monte Carlo+minimization 
algorithm, 44,45 this algorithm includes exchange moves 
between particles at fixed composition as well as ex- 
change moves between configurations at different com- 
positions. Tests on simple Ar n Xei3_„ and Ar„Xei9_ ra 
clusters show that the method is quite efficient, in ad- 
dition to being easy to implement. For these systems, 
we have found that the choice of a very low tempera- 
ture works best as it allows some significant time to be 
spent for optimizing the search for homotops on a same 
common lattice. 

Putative global minima for Ar38_„Xe„ and Kr38_„Xe„ 
clusters have been investigated for all compositions. The 
structure of Ar-Xe compounds is mainly polytetrahedral, 
except at very low doping rates. Kr-Xe clusters not only 
remain as truncated octahedra, but mixing the two rare 
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gases even favors these cubic structures over icosahedra. 
We see some significant trend toward core/surface phase 
separation in Ar-Xe clusters with n > 20 and in all Kr- 
Xe clusters. However, these demixing behaviors are not 
due to the same factors, as Xe atoms favor outer sites 
to reduce strain in Ar-Xe icosahedra, while they occupy 
interior sites to maximize the number of bonds in Kr- 
Xe truncated octahedra. Conversely, Ar38_„Xe„ clusters 
with n < 20 exhibit a higher degree of mixing. Analysing 
the strain in these stable structures confirms the presence 
of a structural transition near n = 20 in these systems. 

Polytetrahedral morphologies were also found as the 
most stable structures of many mixed Ar-Xe clusters 
with 55 atoms, as soon as the relative number of Xe 
atoms was large enough. The general conclusion thus 
seems that the extra strain introduced by mixing these 
different elements penalizes the highly ordered (cubic or 
2-layer icosahedron) structures. 

Within the harmonic superposition approximation, we 
have estimated the temperatures required by the 38-atom 
Kr-Xe clusters to undergo a structural transition toward 
the icosahedral funnel, or toward other octahedral homo- 
tops. For compositions with a doping rate higher than 
3/38, the structural transition temperature was seen to 
occur at temperatures higher than the extrapolated melt- 
ing point. This mainly reflects the special stability of the 
octahedral structures, and has the probable consequence 
that actual melting points increase somewhat. These pre- 
dictions could probably be checked with numerical simu- 
lations. For most compositions, the transitions between 



different homotops of the truncated octahedron are seen 
to be potentially induced by relatively small tempera- 
tures. Therefore particle exchange moves will be neces- 
sary in order that simulations remain close to ergodic. 

Following previous results by other researchers ^2^2^ 
we have found some further evidence that criteria for 
glass formation in bulk materials may also lie in the pa- 
rameters, which are responsible for stable cluster struc- 
tures. Since the atomistic simulation of the dynamical 
vitrification process can generally be much harder than 
obtaining stable configurations of atomic clusters, we ex- 
pect the approach followed in the present theoretical ef- 
fort to be also useful in the community of glasses and 
supercooled liquids. 

The method is obviously not limited to rare-gases, and 
its application to other compounds, especially metallic 
nanoalloys, should be straightforward, except maybe for 
fine tuning its intrinsic parameters. From a methodolog- 
ical point of view, it could also be applied to materials 
with more than two components. Work on ternary sys- 
tems is currently in progress. 
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TABLE V: Global optimization results for Kr38- n Xe n . The energies are given in LJ units for argon, and the symmetry and 
mixing ratios defined by Eq. Q are reported. 
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